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ABSTRACT 

The mesoscale model Meso-NH is used to simulate the optical turbulence at Mt Gra- 
ham (Arizona, US), site of the Large Binocular Telescope. Measurements of the Cff- 
profiles obtained with a Generalized Scidar from 41 nights are used to calibrate and 
quantify the model ability in reconstructing the optical turbulence above the site. The 
measurements are distributed over different periods of the year permitting us to study 
the model performance in different seasons. A statistical analysis of the simulations is 
performed for all the most important astroclimatic parameters: the C^-profiles, the 
seeing e, the isoplanatic angle Oq and the wavefront coherence time tq. 

The model shows a general good ability in reconstructing the morphology of the 
optical turbulence (shape of the vertical distribution of the Cff) as well as the strength 
of all the integrated astroclimatic parameters. The relative error (with respect to 
measurements) of the averaged seeing on the whole atmosphere for the whole sample 
of 41 nights is within 9.0%. The median value of the relative error night by night is 
equal to 18.7% so that the model still maintains very good performances. Comparable 
percentages are observed in partial vertical slabs (free atmosphere and boundary layer) 
and in different seasons (summer and winter) . We prove that the most urgent problem, 
at present, is to increase the ability of the model in reconstructing very weak and 
very strong turbulence conditions in the high atmosphere. This evidence in the model 
mainly affects, at present, the model performances for the isoplanatic angle predictions 
for which the median value of the relative error night by night is equal to 35.1%. No 
major problems are observed for the other astroclimatic parameters. A variant to the 
standard calibration method is tested but we find that it does not provide better 
results confirming the solid base of the standard method. 
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1 INTRODUCTION 

Mt Graham International Observatory is located in south- 
eastern Arizona, USA, at 3200 m above sea level. The obser- 
vatory consists of three telescopes: the Large Binocular Tele- 
scope (LBT - with two 8.4 m mirrors), the Heinrich Hertz 
Sub-millimetre Telescope (SMT - D=10 m) and the Vati- 
can Advanced Technology Telescope (VATT - D=1.83 m). 
Some studies aiming to characterize the optical turbulence 
above Mt Graham have been done in the past with mea- 
surements mainly retrieved from a Gener alized Scidar (GS) 
lEgner et al.ll2007l : [Masciadri et al.ll2010l '). In this paper we 
investigate the possibility to characterize and predict the op- 
tical turbulence (vertical distribution and integrated values) 
at Mt Graham using an atmospheric mesoscale model called 
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Meso-NH. In a previous paper (|Hagelin et al ."201(3^ we have 
used the same model to investigate the possibility to pre- 
dict the vertical wind speed distribution at Mt Graham. We 
proved that the Meso-NH model provides reliable estimates 
of the vertical distribution of the wind speed at all heights 
from the ground up to 20 km. This wind speed can be used 
therefore for the calculation of the wavefront coherence time 
(ro) on the Mt Graham summit. 

The Meso-NH model has already been used to study the 
optical turbulen ce (OT) at different ast ronomical sites. It 
was first used bv lMasciadri et al.l l|l999al lbh who also devel- 
oped the code for the optical turbulence in the Meso-NH, the 
so called Astro-Meso-NH package including the algorithms 
for the C% parametrization. These were, to our knowl- 
edge, the first C|f-profiles ever simulated with a mesoscale 
model in an astronomical context. In those studies the model 
was proved to be sensitive to orographic effects and to be 
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able to reconstruct C^-profiles well correlated to measure- 
ments provided by a Scidar. Also it was able to discrim- 
inate between the worst and the b est seeing of the mea- 
surement campaign. iMasciadri et aLl|l999a .b) discussed the 
main limitations encountered in the C% modelling and pro- 
posed methods to overcome them. 

Some improvements in the model reliability have been 
achieved more recently than ks to the method for th e mode l 
calibration introduced by IMasciadri fc Jabouilld (|200lh . 
Such a calibration reduces some systematic errors and it 
is based on the tuning of a free parameter called minimum 
turbulence kinetic energy (Emin)- The authors proved that 
the C% can be related to Emin and the calibration aims 
to fix the va lue of Emm using as a refere nce the measured 
C^-profiles. IMasciadri fc Jabouilld (|200ll ') proved that the 
calibrated Meso-NH model reduced some of the previous 
systematic errors in the model and the shape of the result- 
ing C^-profiles fitted better with measurements than the 
uncalibrated model. The statis tic reliability of t his m ethod 
has been proved only later by IMasciadri etal] (120041 ) who 
used the calibrated Meso-NH to simulate the optical turbu- 
lence at the San Pedro Martir Observatory on a sample of 
10 nights, for which there were C^-profiles obtained with 
a Generalized Scidar and micro-thermal sensors mounted 
on radiosondes. Measurements provided by different instru- 
ments permitted the authors to prove that the dispersion 
between measurements and simulations was of the same or- 
der as the dispersion between measurements provided by 
different instruments. This qualified the numerical technique 
and the Meso-NH model as potentially useful to perform au- 
tonomous estimates of the optical turbulence. More recently 
(IMasciadri fc Egnedl2006l l the Meso-NH has been used in an 
autonomous way (after calibration) on a sample of 80 nights 
to simulate the optical turbulence at San Pedro Martir. 

The algorithms for the op tical turbulence parametriza- 
tion introduced by IMasciadri et al.l (Il999al) as well as the 
calibration procedure ( Masciadri Sz Jabouilld I2001I ) have 
been late r on implemented i n other mesoscale models su ch 
as WRF (|Klemp et al.ll2008l : iBusinger fc Cherubinill2010l '). 

The goal of this paper is to use the Meso-NH model 
to simulate the optical turbulence at Mt Graham using, 
as a reference, a sample of measurements associated to 
43 nights. A sub-sample of this rich statistical sample is 
used to calibrate the model. The outputs of the model are 
then compa red to measurements obtained with a General- 
ized Scidar (|Masciadri et al.|[2010l ). The sample of measure- 
ments is around four t imes larger than the previous samples 
(|Masciadri et al.ll2004l ') and the largest ever used at present 
for this purpose. For this reason we can: (1) verify if the 
calibration permits the model to still obtain a typical verti- 
cal distribution well correlated with measurements and/or if 
this correlation increases/decreases; (2) study how the cor- 
relation between numerical calculations and measurements 
is deteriorated when we consider in the sample data not used 
for the calibration; (3) study the ability of the Meso-NH 
model in reproducing the seasonal differences in the opti- 
cal turbulence. The measurements are, indeed, quite evenly 
distributed over different periods. In this paper we also test 
a variant of the meth o d of model calibration proposed by 
IMasciadri fc Jabouilld (|200lh . hereafter MJOl. 

In Section [2] we briefly present the measurements used 
as a reference in this study. In Section [3] we present the 



Meso-NH model and describe the model configuration used 
in this study. In Section |4] we describe the calibration pro- 
cedure. In Section [5] we present the results of this study in 
four subsection dedicated respectively to the vertical distri- 
bution of the optical turbulence (C^-profiles), the seeing e, 
the isoplanatic angle 9o and the wavefront coherence time 
tq. In Section |6] the conclusions of this study are presented. 



2 MEASUREMENTS OF REFERENCE 

C% profiles obta i ned with a Generali zed Scidar (GS) 
(jEgner et al.l l2007l : IMasciadri et al.l l2010l l mounted on the 
VATT have been used as a reference for this study. The GS 
is based on the observation of binaries with a typical sepa- 
ration 6 within (3-10) arcsec, the binary magnitude mi,m2 
< 5 mag and A(mi,2) < 1 mag. The GS needs a telescope 
with a pupil size > 1.5 m. The C|^-profiles are obtained from 
the inversio n of the auto-correlation of the scin tillation map 
of binaries (|Fuchs et a]||l998l : lAvila et al.lll997l '). This instru- 
ments provides a vertical distribution of the optical turbu- 
lence on the whole 20 km fro m the groun d with a vertical 
resolution that scales as CTSy^ X{h — hgs) / 9 where h is the 
height from the ground and hgs is the height under ground 
at which the conjugated plane is optically placed. The GS 
provides a vertical resolution that, considering the param- 
eter space just described, is typically of the order of 1km 
on the whole atmosphere, reaching the best resolution near 
the ground (order of some hundreds of meters). GS mea- 
surements are the best choice (with respect to other vertical 
profilers) for a mesoscale model validation for a few rea- 
sons: (1) measurements are obtained with a remote sensing 
technique and they are therefore available for an extended 
period of time during a night, (2) the C% is obtained with a 
completely independent and auto-consistent method requir- 
ing no calibration (3) one can access directly the C% i.e. the 
prime element from which all the integrated astroclimatic 
parameters can be calculated autonomously. Appendix A 
reports the analytical equations that describe how the see- 
ing, the isoplanatic angle and the wavefront coherence time 
are calculated from the C%. 



3 MODEL CONFIGURATION 

The Meso-NH is a non-hydrostatic mesoscale model de- 
veloped jointly by the Centre National des Recherches 
Meteorologiques (CNRM - Meteo-F rance) and Laborat oire 
d'Aerologie in Toulouse (France) (|Lafore et al.l Il998l ). It 
is a grid point model based on the anelastic approx- 
imation to efficiently fi l ter ou t the acoustic waves. A 
iGal-Chen fc SommerviUd (|l975l ) coordinate on the vertical 
and a C-grid in the formulation of Arakaw a et al.l (| 19761 ) for 
the spatial digitalization is used. The temporal scheme is an 
explicit three- time-level leap-frog scheme with a time filter 
(|Asselinlll973 ). The turbulent scheme is a one-dimensional 
1.5 closure scheme (|Cuxart et al. I I2OOOI ). The model per- 
mits the use of different mixing leng t hs. In this paper we 
used the 1-D iBougeault fc Lacarrerd (|l989l ) mixing length 
(BL89). The surface exchanges are computed in an ex- 
ternalized surface scheme (SURFEX) including the physi- 
cal package ISBA (Interactions Soil Biosphere Atmosphere) 
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Figure 1. Topographic areas covered by the models. The left hand panel is the outermost model (Model 1) with a resolution of 10 km. 
The black square indicates the location of the sub-model (Model 2). The central panel is the sub-model with 2.5 km resolution. The 
black square show the location of the innermost model. The right panel shows the innermost model (Model 3) with a 500 m resolution. 
The colour-scale indicates the height above sea level. 



Table 1. Meso-NH configuration. The second column gives the 
horizontal resolution AX of the imbricated models, the third col- 
umn the grid-points, the fourth column the surface covered by the 
models and the fifth column the time steps used by the model. 





AX 


Grid points 


Surface 


Time step 




(km) 




(km) 


(s) 


Model 1 


10 


80 X 80 


800 X 800 


30 


Model 2 


2.5 


64 X 64 


160 X 160 


6 


Model 3 


0.5 


120 X 120 


60 X 60 


3 



l|Noilhan fc Planton|[l989l ) that controls the air/ground tur- 
bulent fluxes budget of Meso-NH. The model can simulate 
the temporal evolution in three dimensions of the classical 
meteorological parameters, such as wind speed and direc- 
tion, potential temperature, p ressure, and so on. Mes o-NH 
uses the code implemented bv lMasciadri et al] (|l999al lb[) to 
forecast the optical turbulence {C% 3D maps) and all the as- 
troclimatic parameters deduced from the 0%. We will refer 
to the 'Astro-Meso-NH code' to indicate this package. The 
integrated astroclimatic parameters are calculated integrat- 
ing the C% with respect to the zenith in the Astro-Meso-NH 
code. 

The Meso-NH model is run, in this paper, in a grid 
nesting mode using three imbricated models with different 
resolutions (Table [l}. They are all centred on the Mt Gra- 
ham Observatory (32.7013°N, 109.8919° W). The size of the 
outermost model (Model 1) is 800x800 km which covers most 
part of south-eastern Arizona, south-west New Mexico and 
also a part of north-western Mexico, see Fig. [T] The middle 
model (Model 2) covers 160x160 km, using a grid-size of 2.5 
km. The innermost model (Model 3) has a resolution of 500 
m and is covering an area of 60x60 km. The vertical reso- 
lution is the same for all three models with 54 vertical grid 
points, reaching up to 20 km above the ground. The first 
vertical grid point is located 20 m above the ground and 
thereafter the grid is determined by a logarithmic stretch- 
ing (20%) up to 3.5 km above the ground. Above 3.5 km the 
resolution is almost constant and equal to ~600 m. 

The model is initialized and forced every 6 hours at the 
synoptic hours (00:00, 06:00, 12:00 UTC) with analyses from 
the the ECMWF (European Centre for Medium-Range 



Weather Forecasts). The model runs for 12 hours, the first 
two hours of every simulation are rejected because the 
model is still adapting to the orography. The output from 
the remaining 10 hours (19:00 - 05:00 Lxfl is used for the 
characterization of the optical turbulence at Mt Graham. 
The Astro-Meso-NH model provides the vertical profile of 
the C% for every two minutes at the grid point located 
at the astronomical Observatory. Fig. [2] (left) shows an 
example of the temporal evolution of the C% extended on 
20 km obtained with the Astro-Meso-NH package. Fig. [5] 
(right) shows the correspondent temporal evolution of 
the me asured C^^-profil e s obt ained with the Generalized 
Scidar (|Masciadri et al.| |2010|) . In these figures we can 
appreciate the characteristics of the simulated C^-profiles. 
Most of the turbulence layers at the different heights of 
the atmosphere appear well reconstructed by the mode{3. 
The spatial variability of the simulated turbulence is in 
general smoother than the measured one. This effect is more 
evident in the high part of the atmosphere. This is what we 
expect from the calculations of a parametrized parameter 
that is not explicitly resolved. Also the temporal variability 
of the simulated turbulence appears smoother with respect 
to the observed one. It is difficult to say that the value of 
a parametrized parameter can be predicted at a precise 
time t=to. For this reason so far we preferred to provide 
averaged estimates of the optical turbulence. We calculate 
the mean of the -profiles simulated of the whole night 
(with a sampling of the C|^-profile every 2 minutes in the 
interval 19:00 - 05:00 lt) and we associate the result to the 
mean of the observed C^r-profiles obtained during the same 
night with the Generalized Scidar. We compare therefore 
typical measured and simulated C^-profiles. Hereafter we 
discuss the performance of the Meso-NH model under this 
assumption. This logic is the same used in all the previous 
studies done on the modelling o f the o p tical turbulence 



with Mes o-NH (iMasciadri et all 1119991: 
(|200lD . IMasciadri fc Garfiad (l200lD 



Masciadri et al 



Masciadri et al 



^ LT = Local Time. 

^ Near the ground the shape of the is necessarily different 
because the Generalized Scidar produces the typical bump due to 
the limited vertical resolution that spread the turbulent energy 
also under-ground. 
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Figure 2. Left: temporal evolution predicted by the Meso-NH model on 1 March 2008 from 00 to 12 UTC. The simulated C^-profiles 
extend on 20 km. The temporal sampling is of one C|^-profile each 2 min. Right: C|^-profiles measured by the Generalized Scidar during 
the same night. 



(|2004). iMasciadri et all (l2004h. iMasc adri fc Egneij (j2006h . 
iLascaux et al.l l|2009l 'l. lLascaux et al.l (|201(]| '1'1. 



4 MODEL CALIBRATION 
4.1 MJOl method 

The method for the model cal ibration implemented in Meso- 
NH has been proposed by IMasciadri fc Jabouillel (l200lh 

i herea fter MJOl) and later on validated by Masciadri et al.l 
2004 ). It is based on the idea that the model depends on 
a free parameter (E^in i.e. the minimum kinetic energy) 
that can be considered as a sort of background climatologi- 
cal noise. In the regions in which the dynamic turbulence is 
well developed the model rapidly forgets the Emi„ value and 
the kinetic energy overcomes this value. The 0% is therefore 
not dependent on the E,„i,i in these regions. However, in 
stable regions, the authors proved that the C% cc E^f^. In 
each C|f-profile it is possible to identify typically more than 
one vertical slab in which turbulence is in stable regimes 
(MJOl-Fig. 1). If we change the value of Emi„ it means that 
we assume that the atmosphere has more or less inertia in 
different vertical slabs and the thermodynamic instabilities 
require more or less energy to trigger turbulence in some 
regions of the atmosphere. As a consequence the threshold 
(the seed Emin) should be different in the respective verti- 
cal slabs. To identify the optimal values of Emi„ the main 
idea in the MJOl method was therefore to divide the whole 
atmosphere (~ 20 km) in a finite number (5-6) of vertical 
slabs and to optimize the value of Emin in each vertical slabs 
minimizing the differences between simulated and measured 
C|^-profiIes i.e. minimizing the function for each night 



2 



JVfc 



\0'm,k ' ^m,i ym,i\ 



(1) 



where Nfc is the number of levels in the vertical slab k, ym,i 
is the average of the measured C% sampled on Nfe levels 
for each night m, Xmj is the simulated C% sampled on 



levels for each night m and am,k is the free coefficient that 
has to b e fixed minimizing the functi on Xm k- We refer the 
reader to lMasciadri fc Jabouillel l|200lf l for details, we simply 
summarize here that, after an average of all the am,k with 
respect to the number of nights m we obtain a coefficient 
for each vertical slab. Knowing the value of ak for k — (1,K) 
(with K the number of vertical slabs), the kinetic energy is 
modified as: 



E* . 



Err 



3/2 



(2) 



Once identified the optimized Emin, characterized by K 
steps, this is implemented in the model that is run again 
for each night. We highlight that the same new Emin is used 
to simulate the new C% for all the nights. This is what we 
call 'an output obtained with a calibrated model' and these 
are the results discussed in the next section, with respect to 
measurements. 

The MJOl method reduces some systematic error and 
reconstructs a typical mean C^-profile better in agreement 
with mea surements. The method has been validated statis- 
tically by IMasciadri et all (|2004h on a sample of 10 nights 
using measurements provided by different instruments and 
taken simultaneously. It is worth to note that the method 
optimizes Emin with respect to the mean C^-profiles simu- 
lated on the total sample used to calibrate the model. The 
calibration is based basically on the conservation of the tur- 
bulent energy (J=C^ ■ AH) in each vertical slab. The dif- 
ferences between measurements and simulations calculated 
in each night, are in general, obviously larger than the dif- 
ference between the statistical values (average in this case) 
as we will discuss later in the paper. The interest in increas- 
ing the statistical sample for the calibration is to verify how 
the reliability of the model changes increasing the number 
of nights. 

In this study the model calibration started with a se- 
lection of the sample of nights on which we calibrated the 
model. A brief digression is necessary here. In an ideal case, 
assuming to have a very rich sample of measurements (for 
example a year of measurements), one should calibrate the 
model using as many as possible different nights in order to 
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approach a profile of Emi„ characteristic of the site. Besides, 
once the model is calibrated, the ideal case should be to have 
an equally rich sample of independent measurements (and 
associated simulations) to investigate the performance of the 
model after calibration. It is also worth to highlight that the 
goal of the calibration is to reduce/eliminate a precise sys- 
tematic error and therefore it is our interest to eliminate 
from the sample used for the calibration, all cases that cor- 
responds to unusual results that might be associated to a 
failure of the model due to whatever reasons. These data, if 
introduced in the calibration sample, might bias the calibra- 
tion. In other words, for the calibration a subjective selection 
is not only allowed but it is suggested without diminishing 
the confidence in the results. 

In this study we can access to a sample of 43 nights. In 
spite of the great number of nights with respect to previous 
studies, it is still a too small number to be able to perform 
the strategy just described using two independent samples 
equally statistically rich. We therefore decided for an alter- 
native solution that is a sort of compromise that maximize 
the progresses (answers to open questions) we can do in this 
field with this data-set. First, we are interested in limiting 
the sources of uncertainties in our sample of reference. The 
assumption done is necessarily to consider the measurements 
as a reference. In other words, we do the approximation that 
measured C^-profiles match with the true C^-profiles. For 
this reason we eliminated definitely from the original sample 
of 43 nights, two nights (March 2 and 3, 2008) because the 
relative shape of the mean of the C^-profile in each night 
is characterized by unusual features and these could affect 
the calibration process. Among the remaining 41 nights, we 
subtracted a sub-sample of seven nights (~ 17% of the total 
sample) in which the model provided unusual results that 
might let us think that for some reason the model did not 
work correctljjfl independently from the algorithm used for 
the C% parametrization. We therefore calibrate the model 
on the remaining sample of 34 nights that is the richest sam- 
ple achievable with this data-set. We study therefore the 
statistical performances of the model calibration done on 
34 nights separating the calibration aspect from the model 
score of success. After that, we quantify the model perfor- 
mances on the whole sample of 41 nights i.e. a sample includ- 
ing also the values of the seven nights that have not been 
considered in the calibration. As said before, this sample 
contains many cases that might be associated to bad model 
performances. This strategy permits us to do some steps 
ahead. We achieve, at the same time: (i) to optimize the 
calibration process, (ii) to discuss the model performance 
including an independent sample of nights not considered 
in the calibration and (iii) to investigate the model perfor- 
mance when it contains a realistic sample of failure cases. 



^ It is not important to discuss the causes for this potential dis- 
crepancies because we do not discard these nights from the total 
sample of nights but simply from the cahbration sample because 
they might bias the calibration process. The typology of "un- 
usual" results are for example, an overestimate of the turbulence 
near the ground. Whatever is the cause for this potential discrep- 
ancy (initialization data not representative of the atmospheric 
flow in a particular night or a specific failure of the model) wo 
simply consider that this is what can happen in a real operational 
prediction. 



After the identification of the samples, we divided the 
atmosphere in six vertical slabs above (0 - 400 m), (400 
m - 2 km), (2 - 7 km), (7 - 11 km), (11 - 13 km), (13 - 
20 km) and we applied the model calibration for h > 400 
m. Below this height, we observed that the model does not 
depend any more on Emin and the seed is quickly forgotten 
by the model during the simulation. As already explained 
in previous papers, the thickness and the number of the 
vertical slabs is arbitrary. We selected regions in which the 
turbulence seems to be characterized by similar trends. 

4.2 Variant to MJOl method 

In this paper we tested if the calibration procedure for Emin 
done on each model level instead of vertical slabs (as is the 
case of the MJOl method) might produce substantially bet- 
ter results or not. The minimization of the function is 
done for each model level: 

E^in,i = Emm ■ aj"^ i = 1,N. (3) 

where N is the total number of the model levels. This 
variant has been suggested by several colleagues in pri- 
vate communications after the pub lication of first r esults 
obtained with the MJOl method l|Masciadri et al.l |2004| : 
iMasciadri fc Egneij[2006 '). From a pure mathematical point 
of view such a method might in theory work better on the 
calibration sample (34 nights) because the mathematical fit 
has more constraints (a total number of measurements equal 
to the number of levels) on which to constrain the fit. How- 
ever, from a physical point of view this method is charac- 
terized by a questionable physical assumption. It is indeed 
hard to justify it because it would be as to admit that the 
climatological noise Emin is different at each model level but 
the values of Emin are the same in each model levels for all 
the nights. This is the reason why we think that, at present, 
the original MJOl version is the most solid approach from a 
physical point of view. However, in this paper, taking advan- 
tage of the rich statistical sample, we decided, independently 
from the arguments that might justify (or not) the MJOl 
variant, to simply test this variant to verify, first, if the gain 
is effective and, second, to check if this gain is maintained 
with samples containing nights not included in the calibra- 
tion sample (in our case the 41 nights sample). To simplify 
the discussion we call hereafter this method MJOl*. 



5 RESULTS 

To discuss the results we divided the solar year in two peri- 
ods: the summer (April - September) and the winter (Octo- 
ber - March). We investigate the optical turbulence vertical 
distribution (C^-profiles) and the three major integrated as- 
troclimatic parameters (seeing, isoplanatic angle and wave- 
front coherence time). For the C% and the seeing we show 
first the results obtained with the sample of 34 nights used 
for the calibration and, afterwards, the results obtained with 
total number of 41 nights. For the isoplanatic angle and the 
wavefront coherence time we treated directly the total sam- 
ple of 41 nights. 

When simulated versus measured C|^-profiles are 
treated it is important to define how to statistically analyse 
the data. This is not univocal. As an example. Fig. [3] shows. 
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Summer Winter 




Figure 4. The average C^-profiles measured by the Generahzed Scidar (red Une) and the Meso-NH model (hght blue line: MJOl; dark 
blue line: MJOl*) after calibration. Left: the average for the total sample of 34 nights. Centre: same for the summer (14 nights). Right: 
same for the winter (20 nights). The green line is plotted in the (400 m - 17 km) range and it shows the average C^-profile obtained by 
the model before the calibration. 
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Figure 5. The average C|^-profiles measured by the Generalized Scidar (red line) and the Meso-NH model (light blue line: MJOl; dark 
blue line: MJOl*) on the total sample of measured nights. Left: the average for the total sample of 41 nights. Centre: same for the summer 
(15 nights). Right: same for the winter (26 nights). 
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Figure 3. Observed C|^-profiles. Red Line: median of all the 
C^-profiles (15956) of the total sample of 41 nights. Green line: 
median of the 41 C^-profiles, each one is the average of the 
of one night. Light blue line: mean of all the observed C|^-profiles 
(15956) of the total sample of 41 nights. Dark blue line: mean of 
the 41 C|^-profiles associated to each night and obtained averag- 
ing the of one night. 

for a pedagogic approach, four C^-profiles calculated from 
the sample of observed C^-profiles: (i) the median of all 
the individual CAr-profiles of the whole sample of 41 nights 




10 " 10 " 10 " 10 " 10"'^ 10" 



CnC"! ) 

Figure 6. Minimum and maximum values of the observed 
and simulated in the sample of the 41 nights. Red lines refer to 
measurements from the GS. Light blue and dark blue lines refer 
respectively to the MJOl and MJOl* models. Values of are 
shown for h > 400 m. 



(15956 profiles) (red line); (ii) the median of the 41 C%- 
profiles associated to each night. Each C|^-profile is obtained 
averaging the C% of a night (green line); (iii) the mean of 
all the individual C^-profiles (15956 profiles) of the whole 
sample of 41 nights (hght blue line); (iv) the average of the 
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Table 2. Calibration sample (34 nights): average seeing in the total atmosphere, boundary layer and free atmosphere. The dome seeing 
is removed from total seeing and the boundary layer seeing. The turbulence contribution provided by the first 20 m from the ground is 
excluded in the contribution of the Meso-NH model. The seeing is given in arcseconds using a wavelength of 500 nm. 



Generalized Scidar Model MJOl Model MJOl* 

£tot £bL £FA £tot SbL £FA £tot £bL ^FA 



Total (34 nights) 0.69 0.50 
Summer (14 nights) 0.46 0.26 
Winter (20 nights) 0.85 0.66 


0.41 
0.34 
0.46 


0.66 0.45 
0.52 0.29 
0.76 0.56 


0.41 
0.37 
0.43 


0.74 
0.60 
0.84 


0.49 
0.32 
0.60 


0.48 
0.47 
0.49 




. Relative error (measured in %) calculated for the calibration sample and 


the total sample. 








Calibration sample: 34 nights 






Total sample: 41 


nights 






Model MJOl 


Model MJOl* 


Model MJOl 




Model MJOl* 


£tot ^BL ^FA ^tot 


£bl 


£FA £tot 


£bl 


sfa 


£tot 


SBL 


£fa 


Total 3.6 9.2 0.4 8.2 
Summer 13.5 11.6 13.0 31.9 
Winter 10.1 14.9 6.3 0.7 


2.1 
21.1 
8.4 


17.6 7.2 
36.8 13.8 
7.6 5.1 


9.0 
10.3 
8.7 


3.2 
13.9 
10.4 


18.0 
32.3 
13.6 


15.4 
19.5 
14.4 


13.2 
37.8 
2.9 



41 C^-profiles associated to each night. Each C|,-profile is 
obtained averaging the Cjf of a night (dark blue line). The 
difference between these C^-profiles tells us that, depend- 
ing on how we treat statistically the data we obtain different 
results. In studies on site characterization done with mea- 
surements taken with monitors approach (i) is frequently 
used. Our context, however, is different. 

As we said, it is hard to predict a C% at a precise time 
(t=t*) of the night and, at present, to investigate quanti- 
tatively the performances of the model our goal is to asso- 
ciate the average of the C|^-profiles simulated during one 
night with the average of the measured C^-profiles of the 
same night. For this reason we are forced to eliminate the 
approach (i) and (iii). Moreover the calibration is based ba- 
sically on the fitting of the average of M C|^-profiles (M 
is the number of nights). We selected therefore the criterion 
(iv) and we considered the 'mean' as a statistical operator to 
quantify the mod el performances o f the C %. This is the same 
approach used in iMasciadri et al.l (|2004l ) . For homogeneity 
we considered the same criterion with the integrated astro- 
climatic parameters (seeing, isoplanatic angle and wavefront 
coherence time). 



5.1 C^-profiles: optical turbulence vertical 
distribution 

Figure [4] shows the average of the C^-profiles measured with 
the Generalized Scidar and simulated with the calibrated 
model with the methods MJOl and MJOl*. In the same pic- 
ture is reported also the average C% before the calibration 
as well as those calculated in the two seasons: summer and 
winter. Figure [5] shows the same for the 41 nights. 

The agreement between the morphology (shape) of the 
averaged C^r-profiles measured by the Generalized Scidar 
(red lines) and calculated with the Meso-NH (light blue 
lines: MJOl; dark blue lines: MJOl*) is very good when look- 
ing at the average of aU the 34 nights (Fig.Hleft). The model 
can reconstruct all the major typical features of the C%- 
profile such as the position (at ~ 10 km) and shape of the 
secondary peak. Between 4 and 8 km above the ground the 
MJOl* method seems slightly better than the MJOl method. 



The vertical distribution of the optical turbulence in the win- 
ter (Fig.U- right) is also very well described by the Meso-NH 
model. However in summer (Fig. [4]- centre) the model seems 
to overestimate the C% in the free atmosphere. However, the 
MJOl method is better than the MJOl* method. The MJOl 
method is indeed well correlated to measurements up to 
around 8 km from the ground while the MJOl* method over- 
estimates the C% starting from 4 km from the ground. The 
'a effect' (|Masciadri fc EgneJlioO^ : IMasciadri et al.ll2010l ) 
is well reconstructed by the model too: in summer the sec- 
ondary peak shifts to a higher altitude as shown by the mea- 
surements. At the same time the strength of the secondary 
peak slightly decreases as observed in measurements. No 
major differences are appreciated in the shape of the av- 
eraged C^-profiles between the calibration case (34 nights) 
and the whole sample (41 nights). Therefore the inclusion in 
the sample of some more nights does not seem to produce 
a major impact on the morphology (shape) of the averaged 
C^r-profile. Which is the cause of the model overestimate 
of the C% in the free atmosphere in summer? Looking at 
Fig. [6]we note that, in the free atmosphere, the model vari- 
ation between the minimum and the maximum values of 
C% is substantially smaller than what observed with mea- 
surements. The higher the altitude, the higher is the model 
inertia and, as a consequence, the model has more difficulties 
in reconstructing the extreme values (minimum and maxi- 
mum) in the high part of the atmosphere. One could expect 
therefore a slight underestimate of the model in winter and 
a slight overestimate in summer. The fact that a model over- 
estimate is observed only in summer can be due to the fact 
that, in this season the number of nights for the 34 and 41 
nights samples are, respectively, 14 and 15 nights while in 
winter we have 20 and 26 nights respectively. It is therefore 
possible that the winter data got a more important weight 
in the calibration. Also in summer we got all nights with 
extremely weak turbulence in the free atmosphere. With a 
richer and more homogeneous statistical sample we should 
therefore expect a decreasing in strength of the discrepancy 
and a more symmetric discrepancy with respect to the two 
seasons. We do not think that there is a specific problem 
for the summer period. Nevertheless it is a fact that the 
model variability needs to be improved in the high part of 
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Table 4. Total sample (41 nights): average seeing in the total atmosphere, boundary layer and free atmosphere. Same as Table [2] 

Generalized Scidar Model MJOl Model MJOl* 

etot £bl ^FA ^tot £bl ^FA £tot sbl ^fa 

Total (41 nights) 0.71 0.51 0.43 0.76 0.56 0.41 0.84 0.59 0.48 

Summer (15 nights) 0.45 0.26 0.34 0.51 0.28 0.38 0.60 0.31 0.46 

Winter (26 nights) 0.87 0.67 0.48 0.91 0.73 0.43 0.98 0.76 0.49 



the atmosphere to be able to well detect all the cases of very 
strong as well as very weak turbulence conditions. We are 
working at present on this topic. 



5.2 Seeing: e 

The seeing depends on the C% as Eq IA2I in Appendix A. 
Figure [7| shows the total seeing (etot), the boundary layer 
(ebl) and the free atmosphere (efa) simulated by the model 
(MJOl and MJOl*) plotted against the respective seeing ob- 
served by the GS for the sample of 34 nights after calibra- 
tion. The boundary layer is defined as the first kilometre 
above the surface and the free atmosphere from the top of 
the boundary layer up to 20 km above the surface. All values 
of the seeing are calculated opportunely by subtracting in 
the measurements from the GS the dome seeing contribu- 
tion and by subtracting in simulations the contribution of 
first 20 m of the C^-profiles (the equivalent of the telescope 
height). In Figure [7] (bottom-right of each panel) is reported 
the correlation coefficient calculated for the total sample of 
41 nights. Table [2] and Table [Sj-left report the values of the 
average seeing and the relative error for the calibration sam- 
ple (34 nights). Tableland Table [31-right reports the same 
for the total sample of 41 nights. The correlation coefficients 
are indicated for the MJOl and MJOl* cases in each panel of 
Fig. \7\ The total seeing and the seeing in the boundary layer 
show a good correlation (c.c. = 0.78-0.82) with measure- 
ments in the calibration sample with no major differences 
between the MJOl and MJOl* cases. It is very encouraging 
that the turbulence near the ground, that represents most 
of the turbulence developed in the atmosphere, is well pre- 
dicted by the model. The correlation decreases for the seeing 
in the free atmosphere. The explanation for such an effect 
is that, in the free atmosphere, visibly the model inertia is 
still too high and the parameters predicted by the model 
varies in a smaller range than what the measurements show 
as discussed in Section [5. II 

The model has at present some problems in identify- 
ing the best (the minimum observed efa) and the worst 
(the maximum observed efa) conditions in this part of the 
atmosphere. Looking at Fig. [Jit appears clear that the prob- 
lem concerns the C% at h > 10 km. However, we note that 
the model (with both methods, MJOl and MJOl*) is able 
to reconstruct a weaker efa in summer than in winter as 
we expect and also as is observed (Table [2|. On the sample 
of calibration no substantial and systematic differences are 
noted on results obtained with the MJOl and MJOl* meth- 
ods. The very encouraging result is that the relative error 
for the seeing in all the three regions of the atmosphere 
(Table O-left) is very good even in the free atmosphere. Re- 
sults are equally very good if we look at the individual sea- 
sons. The important conclusion is that the relative error for 



the seeing on the whole sample is within 9.2% for the best 
method (MJOl) and within 17.6% if we consider also MJOl* 
methocQ. If we consider the two sub-samples (summer and 
winter) the relative error decreases slightly still maintains 
good performances. In conclusion we observe relative errors 
remarkably good in spite of a modest correlation in the free 
atmosphere. We also note that the correlation is very good 
for seeing smaller than 1 arcsec and deviations are more 
typical when there is strong turbulence. In conclusion we 
can say that the method MJOl confirms to be effective and 
it improves the model reliability reducing some systematic 
effects. 

If we look at the total sample of 41 nights (Tableland 
Table [31-right), including nights not used for the calibration 
we observe that the correlation decreases for etot and ebl 
(from 0.82 to 0.61). For efa, the correlation decreases in a 
less consistent way and the MJOl method seems better than 
the MJOl* method in this region. Most of the new nights 
included in this sample belong to the winter period and the 
reason for the decrease of the correlation is that the model 
reconstructs a seeing in the boundary layer that is statisti- 
cally too strong. In this case, the MJOl method is visibly bet- 
ter than the MJOl* method in both cases: the total sample 
and the summer and winter periods. Method MJOl achieves 
a remarkably good result with a relative error within 9.0% 
for the total sample and a maximum relative error of 13.9% if 
we look at the sub-samples (summer and winter). The MJOl 
method is also better than MJOl* method in reconstructing 
the weakest and the strongest values in the two season in 
the partial efa and ebl- Also in this case the model recon- 
structs a Efa that varies in a smaller range with respect to 
what observed. Equally to the calibration sample, we can 
however observe that the model can reconstruct a weaker 
seeing in summer than in winter as expected. 

Finally, to study the performance of the model night by 
night in a compact way, we have calculated the cumulative 
distribution of the relative errors with respect to measure- 
ments for the Etot, Ebl and efa obtained with the MJOl 
method and the MJOl* method (Fig. [H]). We note that the 
relative error in the three regions of the atmosphere is abso- 
lutely remarkably good. The median values (50% of times) 
of the relative error is within 26.3% in the three regions of 
the atmosphere for the MJOl method, which is visibly better 
than the MJOl* method in this case. This means that, even 
if the correlation for the efa is not as good as for Etot and 
Ebl, the relative error night by night is maintained within 
21%. The biggest relative error is obtained in the boundary 



* This indicates that increasing the number of night s the cali- 
bration attain more accurate results. iMasciadri et ah I l l2004h on 
a sample of 10 nights obtained a relative error on the average of 
the order of 20%. 
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Figure 7. Values of the seeing from the simulations (MJOl xs, MJOl* circles) plotted against the measurements in the calibration sample 
(34 nights). The left hand panel shows the total seeing, the middle panel the seeing in the boundary layer, the right hand panel the free 
atmosphere seeing. The dome seeing has been removed from the Generalized Scidar data and the seeing in the first 20 m (correspondent 
to the height of the telescope) has been removed from the Meso-NH simulations. In bottom-left corner of each panel the correlation 
coefficient calculated for MJOl and MJOl* for the calibration sample is reported. In bottom-right corner (red) the c.c. calculated for the 
total sample of 41 nights. 




Figure 8. Cumulative distribution of the relative errors of each single night of the total seeing (£tot)i the seeing in the boundary layer 
(sbl) a-iid the seeing in the free atmosphere (efa) for the total sample of 41 nights. Bold line: MJOl method. Thin line: MJOl* method. 



layer. Even if this is the region in which the correlation is the 
best one, the impact of an error of the model can be more 
important than an error produced in the free atmosphere. 

In conclusion we can state that, on the total sample 
of nights, the MJOl method provides globally better results 
than the MJOl* method. This means that the MJOl method, 
which is based on a minor number of constraints than the 
MJOl* method, has a better performance probably because 
it is based on more robust physical assumptions. A larger 
number of constraints can be effective on the calibration 
sample but, when an independent sample of nights is taken 
into account, this theoretical advantage seems to lose its 
effect. 

5.3 Isoplanatic angle: 9o 

Figure [S3] shows the simulated versus the measured isopla- 
natic angle 9o for the sample of 41 nights. 6o depends on the 
C% as Ea lA3l in Appendix A. The values reconstructed by 
the model span a much smaller range (1.7-2.2) arcsec than 
what has been observed (0.9-5.6) arcsec. The average of the 
values obtained for all the nights is in good agreement with 



measurements. The average of the simulations is equal to 
2.06 arcsec for MJOl (1.82 arcsec for MJOl*) versus 2.65 
arcsec measured by the GS. The difference between simu- 
lations and measurements is due to a small off-set in the 
model calibration. As we have already discussed in the pre- 
vious sections, the calibration sample is slightly biased in the 
high part of the atmosphere by the fact that measurements 
are not completely homogeneously distributed in the two 
seasons and in the summer we have many nights with ex- 
tremely weak turbulence in the high part of the atmosphere. 
Looking at Fig. [5]-left we can appreciate a generally very 
good agreement between measurements and simulations but 
a very weak overestimate by the model is visible in the C% 
above 10 km. Oq is very sensitive even to small variation 
of the C% in this region since it scales with C%{h) ■ h^^^ . 
Therefore a very weak difference in the €% in the high part 
of the atmosphere can produce sub-arcsecond differences in 
the isoplanatic angle. The problem of the small off-set in 
the average should disappear if one uses for the model cal- 
ibration a very rich and homogeneous statistical sample of 
nights (one year for example). However, the failure of the 
model in reconstructing the larger variability of 6o quanti- 
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Figure 9. Isoplanatic angle (0o) from the simulations (MJOl 
stars, MJOl* circles) plotted against the measurements related 
to the sample of 41 nights. 




Figure 10. Cumulative distribution of the relative errors of each 
single night of the isoplanatic angle (do) for the total sample of 
41 nights. Bold line; MJOl method. Thin Une: MJOl* method. 



fied by measurements is not due to the calibration but to 
the fact that, at present, the model shows a larger inertia in 
the high part of the atmosphere as discussed in Section TS. II 
This produces a more modest correlation between simula- 
tions and measurements for 9o. As already said we are at 
present working on this issue to improve the model variabil- 
ity. 

Figure \5A\ shows the cumulative distribution of the rela- 
tive error (night by night) for 9o obtained with both methods 
(MJOl and MJOl*). The median value of relative error for 6o 
is slightly larger (35.1 % for the best method, MJOl) than 
what is observed for the seeing described in the previous 
section. 



5.4 Wavefront coherence time: ro 

To depends on the C% and the wind speed vertical pro- 
files as Eq |A5l in Appendix A. How do we calculate tq for 
mecisurements and simulations? For th e wind speed pro- 
files we have showed in previous studies (jEgner et al.ll2007l : 
iMasciadri et al.ll2010l ') that the best solution is to consider 
a composite profile: the wind speed profiles retrieved from 
the ECMWF analyses (i.e. analyses from General Circula- 
tion Model) in the nearest grid point (~12 km away) for h 
above 1 km and the wind speed as retrieved from the GS 



for h below 1 km. This is necessary because the analyses 
are not sampled with enough high horizontal resolution and 
they are not reliable near the gro und. However, in a more 
recent study (jHagelin et al.|[2010l ) it has been proved that 
the wind speed profiles provided by the Meso-NH model at 
the summit of Mt Graham are very well correlated to the 
wind speed estimated by the ECMWF analyse^ above 1 
km and are well correlated to measurements taken with a 
GS and an anemometer near the ground. To calculate the 
To we considered therefore the wind speed reconstructed by 
Meso-NH above the summit. Due to the fact that for both, 
measured and simulated ro, we used the wind speed pro- 
files as retrieved from the Meso-NH, when comparing the 
To we are basically comparing the effects of the simulated 
and measured C% on the simulated and measured tq. We 
consider the average of the wind speed profiles during each 
night and we calculate the ro for each night. Figure [TT]- left 
shows the simulated versus the measured wavefront coher- 
ence time (ro) for the sample of 41 nights. The data-set is 
well distributed along a straight line showing a very good 
correlation (0.95 for MJOl and 0.96 for MJOl*). Up to 10 
msec the values are very well correlated. A very small bias 
is evident for the extremely good values of ro (> 10 msec). 

To verify which part of the C% (low or high atmosphere) 
mainly affects ro in Fig. llll -left. we calculated a ro in the par- 
tial 0-10 km and 10-20 km ranges (Fig. llH -centre and Fig. 1 111 - 
right). We observe that while the simulated and measured 
values are well correlated in the 0-10 km range, the model 
underestimates the ro in the 10-20 km range. However the 
contribution coming from the 0-10 km region has a much 
more important affect on the calculation of ro on the whole 
atmosphere than the contribution coming from the 10-20 km 
region. This does that ro on the whole 0-20 km is very well 
reconstructed by the model globally. This is a confirmation 
that, also for this parameter, the small bias of the 0% in 
the high part of the atmosphere produces some effect on ro. 
However, differently from Oq, this effect is almost negligible 
(Fig.[TT}left). 

To study the performance of the model night by night. 
Fig. l5.4l shows the cumulative distribution of the relative er- 
ror for ro, simulated with both the MJOl and MJOl* meth- 
ods, and calculated, taking the measurements as a reference. 
The median of the relative error is very good for ro, as good 
as for the seeing, 22.5 % for the MJOl method and 21.8 % 
for the MJOl* method. 



6 CONCLUSIONS 

In this study we discussed the abilities of the Meso-NH 
model in simulating the optical turbulence above Mt Gra- 
ham, site of the Large Binocular Telescope. Simulated C%- 
profiles are compared to a large statistical sample of mea- 
sured C|^-profiles related to 41 nights obtained with a Gen- 
eralized Scidar. This large sample of measurements allowed 



^ The authors also showed that the wind speed is uniform on a 
horizontal scale of some tens of kilometres for h > 1 km. The 
wind speed from analyses close to Mt Graham are in agreement 
with radiosoundings launched from Tucson International airport 
(~120 km from Mt Graham) and this guarantees the reliability 
of analyses. 
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Figure 11. Wavcfront coherence time {to) from the simulations (MJOl stars, MJOl* circles) plotted against the measurements related 
to the sample of 41 nights. 




Rel. error {%) 

Figure 12. Cumulative distribution of the relative errors of each 
single night of the wavefront coherence time (to) for the 41 nights. 
Bold hne: MJOl method. Thin hne: MJOl* method. 



us to study the performances of the model cahbration in 
great detail and to investigate the performance of the model 
in different seasons discussing statistically how^ well the 
model reconstruct the C%, the seeing in different regions 
of the atmosphere, the isoplanatic angle and the wavefront 
coherence time. 

Two different methods of model calibrations have been 
investigated. The calibration is done on a sample of 34 nights 
and the model performances are discussed on the whole sam- 
ple of 41 nights. 

The most important results obtained are: 

(i) Wc have proved that the model calibration definitely im- 
proves the model performance. The morphology of the verti- 
cal distribution of the optical turbulence (average of all Cff 
of the 34 nights) matches very well with the corresponding 
measured C^-profile. If we look at the model behaviour in 
different seasons, wc show evidence of a model overestimate 
of C% in the high part of the atmosphere (h > 8 km) in sum- 
mer. This bias is highly probable because the sample of the 
investigated nights is not completely uniformly distributed 
between summer and winter. With a reasonably richer sam- 
ple (of the order of 1 yr) this small bias can be corrected. 
This evidence proves the necessity of a very rich statistical 
sample to carry out an efficient model calibration. 

(ii) We observed that the model basically always recon- 



structs the most important features of the shape of the mea- 
sured C|,-profile but it still has some difficulties in recon- 
structing the very extreme values (very good and very bad 
turbulence conditions) in the free atmosphere. We are work- 
ing at present on improving these model performances. 

(iii) For the model calibration, the total and boundary layer 
seeing reconstructed by the Mcso-NH model are well cor- 
related with measured values with a correlation coefficient 
(c.c.) of the order of 0.78-0.82. The seeing in the free at- 
mosphere is more weakly correlated to measurements (c.c. 
~0.58). The relative errors are however extremely good in 
all the three regions (total seeing, in the boundary layer and 
free atmosphere). The best method MJOl provide a relative 
error within 9.2% on the total sample and within 14.9% if 
we consider the two seasons. These percentages remain basi- 
cally the same if we consider the whole sample of 41 nights. 

(iv) When wc consider the total sample of 41 nights the cor- 
relation between simulations and measurements decreases 
slightly for etot and sbl (c.c. ~0.60), much less for epA (c.c. 
-0.47). 

(v) If we consider the cumulative distribution of the relative 
errors night by night (i.e. the typical conditions of the op- 
erational mode) we find extremely encouraging results. The 
median value of the relative errors is indeed extremely small 
with some small variations for the seeing in the three regions 
of the atmosphere but typical of the order of 20%. Therefore 
even if we look at the most difficult conditions typical of the 
operational mode, the model maintains a reasonable good 
result. 

(vi) For the isoplanatic angle the relative error of the aver- 
age is still very good (— 18%) but the median of the relative 
errors calculated night by night is 35.1%. It is the parameter 
with the poorest performances between s, 6o and to- This is 
not surprising because it is the parameter that, more than 
others, is very sensitive to the turbulence in the high part 
of the atmosphere. 

(vii) The wavefront coherence time reconstructed by the 
model shows a very good correlation with measurements 
(c.c. ~ 0.95). We proved that the small bias in the high 
atmosphere on the C^r produces a negligible effect on tq. 
The first 0-10 km represent the most important contribu- 
tions to the final To value. The cumulative distribution of 
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the relative error calculated night by night is still very good 
with a median value equal to 22.5%. 

(viii) The MJOl method to calibrate the model confirms to 
be a solid method and it shows a general better performance 
of the model than the variant MJOl*, which consists in fit- 
ting measurements and simulations on each model level in- 
stead of vertical slabs of a few kilometers. It is worth to note 
that the calibration of the high atmosphere is more delicate 
because a small discrepancy from the measured C^-profile 
can produce a large discrepancy on some astroclimatic pa- 
rameters (such as 60). 

The general conclusions are that the Meso-NH model is able 
to describe the optical turbulence distribution above Mt 
Graham showing very good performances with small relative 
errors with respect to measurements for most of the astrocli- 
matic parameters. The most urgent problem to be solved is 
the improvement of the model ability in reconstructing the 
very weak and very strong turbulent conditions in the high 
atmosphere, i.e. the ability in reconstructing isoplanatic an- 
gles in better agreement with measurements. Finally we have 
put in evidence the importance to have a very rich statistical 
sample of C|f-profiles to efficiently calibrate the model. In 
the future we would like to access an even much richer sam- 
ples of measurements done preferably with a GS in order to 
be able to consider a completely independent sample for the 
model validation. In the future it would also be interesting to 
test the model performance when it is initialized with fore- 
casts from the ECMWF instead of the operational analyses 
in order to use the Meso-NH model as a part of a system to 
forecast the optical turbulence at Mt Graham allowing the 
implementation of a flexible-scheduling management of the 
LBT. 



APPENDIX A: INTEGRATED 
ASTROCLIMATIC PARAMETERS AS A 
FUNCTION OF THE C^. 

The seeing (e), the isoplanatic angle (So), the wavefront co- 
herence time (ro) are defined as: 
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